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A FORMULA FOR THE DUBROVNIK POLYNOMIAL OF 
RATIONAL KNOTS 

CARMEN CAPRAU AND KATHERINE URABE 


Abstract. We provide a formula for the Dubrovnik polynomial of a rational 
knot in terms of the entries of the tuple associated with a braid-form diagram of 
the knot. Our calculations can be easily carried out using a computer algebra 
system. 


1. Introduction 

Rational knots and links are the simplest class of alternating links of one or two 
unknotted components. All knots and links up to ten crossings are either rational 
or are obtained by inserting rational tangles into a small number of planar graphs 
(see ED- Other names for rational knots are 2-bridge knots and 4-plats. The names 
rational knot and rational link were coined by John Conway who defined them as 
numerator closures of rational tangles, which form a basis for their classification. A 
rational tangle is the result of consecutive twists on neighboring endpoints of two 
trivial arcs. Rational knots and rational tangles have proved useful in the study of 
DNA recombination. 

Throughout this paper, we refer to knots and links using the generic term ‘knots’. 
In mi, Lu and Zhong provided and algorithm to compute the 2-variable Kauffman 
polynomial [6j of unoriented rational knots using Kauffman skein theory and linear 
algebra techniques. On the other hand, Duzhin and Shkolnikov (3] gave a formula 
for the HOMFLY-PT polynomial 5, l3J of oriented rational knots in terms of a 
continued fraction for the rational number that represents the given knot. 

A rational knot admits a diagram in braid form with n sections of twists, from 
which we can associate an n-tuple to the given diagram. Using the properties of 
braid-form diagrams of rational knots and inspired by the approach in [3] (namely, 
deriving a reduction formula and associating to it a computational rooted tree), in 
this paper we provide a closed-form expression for the 2-variable Kauffman polyno¬ 
mial of a rational knot in terms of the entries in the n-tuple representing a braid- 
form diagram of the knot. We will work with the Dubrovnik version of the 2-variable 
Kauffman polynomial, called the Dubrovnik polynomial. Due to the nature of the 
skein relation defining the Dubrovnik polynomial (or the Kauffman polynomial, for 
that matter), deriving the desired closed-form for this polynomial of rational knots 
is more challenging than for the case of the HOMFLY-PT polynomial. 

The paper is organized as follows: In Section[2]we briefly review some properties 
about rational tangles and rational knots, which are needed for the purpose of this 
paper. In Section[3]we look at the Dubrovnik polynomial of a rational knot diagram 
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in braid-form and write it in terms of the polynomials associated with diagrams that 
are still in braid-form but which contain fewer twists. Our key reduction formulas 
are derived in Section[4]and used in Section[5]to obtain a closed-form expression that 
computes the Dubrovnik polynomial of a standard braid-form diagram of a rational 
knot in terms of the entries of the n-tuple associated with the given diagram. We 
finish with an appendix containing the Mathematica® code, written by the second 
named author, that computes the Dubrovnik polynomial of a rational knot diagram 
in braid-form, based on the formulas obtained in Section [3] 

The paper grew out of the second named author’s master’s thesis at California 
State University, Fresno. 


2. Rational knots and tangles 


Rational tangles are a special type of 2-tangles that are obtained by applying a 
finite number of consecutive twists of neighboring endpoints starting from the two 
unknotted arcs [0] or [oo] (called the trivial 2-tangles ) depicted in Figure [lj An 
example of a rational tangle diagram is shown in Figure [2] 


oo 



[ 0 ] 


Figure 1. The trivial 2-tangles [0] and [oo] 



Figure 2. The rational tangle T(4, 3,2) in standard form 

Taking the numerator closure, N(T ), or the denominator closure, D(T), of a 
2-tangle T (as shown in Figure |3| results in a knot or a link. In fact, every knot 
or link can arise as the numerator closure of some 2-tangle (see 0). However, 
numerator (or denominator) closures of different rational tangles may result in the 
same knot. 

Numerator and denominator closures of rational tangles give rise to rational 
knots. These are alternating knots with one or two components. It is an interesting 
fact that all knots and links up to ten crossings are either rational knots or are 
obtained from rational knots by inserting rational tangles into simple planar graphs. 
For readings on rational knots and rational tangles we refer the reader to [DI21[5] 

Emmmmm- 

Conway [2] associated to a rational tangle diagram T a unique, reduced rational 
number (or infinity), F(T), called the fraction of the tangle, and showed that two ra¬ 
tional tangles are equivalent if and only if they have the same fraction. Specifically, 
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D(T) 



N(T) 


Figure 3. The denominator and numerator closures of a 2-tangle T 

for a rational tangle in standard form, T(b\ , b 2 , ■ • •, b n ), its fraction is calculated by 
the continued fraction 

F(T) = [h,b 2 , := h -+- - - 

b 2 + ... + - - 

bn -1 + T~ 

On 

where b\ £ Z and b 2 ,...,b n £ Z\{0}. Proofs of this statement can be found 

br [II [5,'51 !T2]. 

A rational knot admits a diagram in braid form, as explained in Figures [4] and [5] 
We denote by D[b\, b 2 ,..., b n \ a standard braid-form diagram (or shortly, standard 
diagram ) of a rational knot, where bfs are integers. 



Figure 4. Standard braid-form diagram, n odd 



Figure 5. Standard braid-form diagram, n even 

A standard diagram of a rational knot is obtained by taking a special closure of 
a 4-strand braid with n sections of twists, where the number of half-twists in each 
section is denoted by the integer b t and the sign of bi is defined as follows: if i is 
odd, then the left twist (Figure [6]) is positive, and if i is even, then the right twist 
is positive (equivalently, the left twist is negative for i even). In Figures [4] and [5] all 
integers bi are positive. 

Note that the special closure for n even is the denominator closure of a rational 
tangle and for n odd is the numerator closure. 

If a rational tangle T(b\, b 2 ,..., b n ) has its fraction a rational number other than 
0 or oo, then we can always find bi such that the signs of all the bi are the same 
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right twist left twist 

Figure 6. 


(see pang). Hence, we assume that a standard diagram of a rational knot is an 
alternating diagram. 

In addition, we can always assume that for a standard diagram D\b\ 1 b 2 ,. ■ ., b n ] 
of a rational knot, n is odd. This follows from the following properties of a continued 
fraction expansion [61, b 2 ,..., b n ] for a rational number | such that bi > 0: 

• If n is even and b n > 1, then 

[bi,b 2 ,...,b n \ = [b 1 ,b 2 ,...,b n - 1,1] 

• If n is even and b n = 1, then 

[b\,b 2 ,... ,b n ] = [&i, b 2 , ..., b n -1 + 1] 

The statement below is well-known (see for example m- 

Lemma 1. The following statements hold: 

(1) D[-bi, -b 2 ,..., -b 2k + i] is the mirror image of D[b 1: b 2 ,.. ■ Mk+ 1]- 

(2) D[bi,b 2 ,..., b 2k+1 ] is ambient isotopic to D[b 2k+1 , ...,a 2 ,a 1 ]. 

(3) D[bi,b 2 ,, b 2k +i] is ambient isotopic to N(T(bi,b 2 ,b 2k +i))- 

We denote by K{p/q) a rational knot with standard diagram D[bi, b 2 ,..., b 2k+ i], 
with bi 0, where 

P 

- = [h,b 2 ,.. .,b 2k+ 1], gcd (p,q) = 1 and p > 0. 

If all bi are positive then q > 0, and if all bi are negative then q < 0. The integer p 
is odd for a knot and even for a two-component link. It is known that two rational 
knots K(p/q) and K(p'/q') are equivalent if and only \ip = p' and q' = q ±1 (mod p) 
(see PE [ESIGil). 


3. The Dubrovnik polynomial of rational knots 


In [6], Kauffman constructed a 2-variable Laurent polynomial which is an in¬ 
variant of regular isotopy for unoriented knots. In this paper we work with the 
Dubrovnik version of Kauffman’s polynomial, called the Dubrovnik polynomial. 

The Dubrovnik polynomial of a knot K , denoted by P(K) := P(K)(z,a), is 
uniquely determined by the following axioms: 


1. P(K) = P(K') if K and K' are regular isotopic knots. 

\ / I D I \ / 


2. P 
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The diagrams in both sides of the second and third axioms above represent larger 
knot diagrams that are identical, except near a point where they differ as shown. 
We will use the following form for the second axiom; 

P (x)= p (x) + 2 P ()C)- 2 P (X). 

and refer to it as the Dubrovnik skein relation. Moreover, we will refer to the 
resulting three knot diagrams in the right-hand side of the Dubrovnik skein relation 
as the switched-crossing state, the A-state, and the B-state, respectively, of the 
given knot diagram. 

The following statement is well-known and follows easily from the definition of 
the polynomial invariant P. 

Lemma 2. If K is the mirror image of K 1 then P(K)(z,a ) = P(K )(— z, a -1 ). 

For more details about the Kauffman polynomial and the Dubrovnik version of 
it we refer the reader to mm- 

The goal of the paper is to give an algorithm which computes the Dubrovnik 
polynomial of a standard diagram D\b\, 6 2 , ■ • ■, b n ] for a rational knot. We will use 
the following notation: 

V[h,b2, • • ■, 6 „] := P{D[b!, b 2 , ■■■, b n ]). 

We will focus on the case with positive integers 6 ,. The case with negative inte¬ 
gers bi follows from Lemma[l]and Lemma[2j For a standard diagram D[bi , 6 2 , ■ ■ • > b n ] 
of a rational knot, we call the integer n the length of the diagram. 

We consider first a standard diagram of length three, £>[ 61 , 62 , 63 ], where bi > 3 
for 1 < i < 3. We obtain the tree diagram depicted in Figure [7J whose edges are 
labeled by the weights of the polynomial evaluations of the resulting knot diagrams 
obtained by applying the Dubrovnik skein relation at the leftmost crossing in the 
original diagram. 



Figure 7. A tree diagram for £>[ 61 , 62 , 63 ] 


Note that the middle leaf of the tree in Figure [7] is obtained after applying a 
type II Reidemeister move. Moreover, the diagram at the bottom of the left-hand 
branch of the tree can be modified to a standard braid-form diagram, as exemplified 
in Figure[ 8 j We obtain the following recursive relation for V\b\, 62 , 63 ], where bi > 3 
for 1 < i < 3: 

'P[b 1 ,b 2 ,b 3 \ = P\b x - 2 , 62 , 63 ] - za~ bl+1 V[l, 6 2 - 1, 63 ] + z'P[b 1 - 1,6 2 ,6 3 ]. 
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b 2 



62 — 1 



Figure 8 . Obtaining a standard braid-form diagram 


These types of recursive formulas form the foundation of our later work. In general, 
as in the case above, the switched-crossing state reduces a section by 2 half-twists, 
the A-state reduces a section by one half-twist, and the B-state reduces the diagram 
by one section of twists while contributing some power of a. Because of this, the 
remainder of our cases will deal with standard braid-form diagrams that have 1 or 
2 half-twists in the first section of twists (corresponding to 61 ), since these relations 
will completely reduce the first section. In addition, the move shown in Figure [ 8 ] 
will be commonly referred to as the sliding move and not shown in detail for the 
other cases. 

Consider now a standard diagram D[ 1,1, b 3 ] and the tree given in Figure [9j 



Figure 9. A tree diagram for D[ 1,1, 63] 


Observe that we applied a type II Reidemeister move to arrive at the diagram 
representing the middle leaf of the tree. We easily see that the Dubrovnik polyno¬ 
mial for the diagram Z?[l, 1, 63 ] satisfies the following relation: 

V[l, 1, 63 ] = a " 63 - zV[b 3 + 1] + zaV[b 3 \. 


Next we consider the standard braid-form diagram D[2, b 2 ,b 3 \. Using the tree 
for the diagram D[ 2, b 2 , b 3 ] depicted in Figure 10 and the sliding move for the dia¬ 
gram representing the leftmost leaf in the tree, we obtain the following polynomial 
expression: 

V[2, b 2l b 3 ] = a b2 V[b 3 ] - za-'Vll, b 2 - 1, b 3 } + zV[ 1, 62 , 63 ]- 


Lastly, consider the tree for the standard braid-form diagram £>[1, 6 2 , 63 ] given 
in Figure El By applying the sliding move to the first two leaves in the tree, the 
Dubrovnik polynomial of D[l,b 2 ,b 3 \ satisfies the expression given below: 

V[l, b 2 , b 3 ] = V[l, b 2 - 2, 63 ] - zV[ 1, b 2 - 1, 63 ] + za b *V[b 3 }). 


There are a few more cases that need to be considered, namely when b\ or b 2 
are equal to 1 or 2. These cases are treated in a similar way as above. We collect 
these cases and those shown above in the following statement. 
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Figure 10. A tree diagram for D[ 2, b 2 , 63] 



Lemma 3 . The Dubrovnik polynomial of a standard braid-form diagram of length 
3 with positive twists satisfies the following relations: 

i. V[b 1 ,b 2 ,b 3 ) = V[bi - 2 ,b 2 ,b 3 }-za 1 ~ bl 'P[l,b 2 - l,b 3 ) + z'P[b 1 - 1,62,63], for 

b\ ^ 3 ,62 2. 

ii. V[ 2 , 62,63] = a b 2 V[b 3 \ - za~ 1 V[l,b 2 - 1 ,63] + ^[1,62,63], for b 2 > 2 . 

iii. V[l, 62, 63] = V[l, b 2 - 2, 63] - zV[ 1 , b 2 - 1 , 63] + za b 2 V[b 3 ], for b 2 > 3 . 

iv. V[b 3 , 1,63] = V[bi - 2 , l,b 3 ]-za 1 ~ bl 'P[b 3 + l]+zV[b 1 - 1 , 1 , 63], for b 3 > 3 . 

v. ^[2,1,63] =aV[b 3 }- za~ 1 V[b 3 + l} + zV[l, 1 ,63]. 

vi. V[l, 2 ,b 3 ]=V[b 3 + l]-zV[l,l,b 3 ] + za 2 V[b 3 \. 

vii. V[l, 1 , 63] = a~ bs - ^[63 + 1 ] + zaV[b 3 \. 


Note that the algorithm that allowed us to arrive at Lemma [3] involved only the 
two leftmost groups of twists in a diagram of length 3. Therefore, the above cases 
can be generalized to a standard braid-form diagram of any odd length n with 
positive half-twists 6 ^. 

The one extra case we need to consider before we generalize the statement in 
Lemma[3]is the standard diagram of length 5, D[ 1,1, 63 ,1, 65 ], shown in Figure 12 



Figure 12. Standard diagram D[ 1,1, 63, 1, 65] 
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The Dubrovnik skein relation applied to the leftmost crossing in the diagram 
D[l, 1 , b 3 , 1 , 65] yields the following expression: 

V[l, 1 ,63, 1 ,65] = a~ b 3 V[ 1 + 65] - zV[ 63 + 1 , 1 ,65] + zaP[b 3 , 1 , 65]. 

Theorem 1 . The Dubrovnik polynomial of a standard braid-form diagram D[bi, b 2 ,... 
with n odd and positive twists satisfies the following relations: 


V[bi,b 2 ,. 

• • ? bn] 

— P[b\ — 2, b 2 , ■ ■ ■, b 

n] -za^^Vi 1,62 - 1,...,6„] 



+zV[bi — 1 ,62,.. 

• ,6„], for bi > 3,62 > 2. 

V[ 2 ,b 2 ,- 

■ • ? bn] 

= a b 2 V[b 3 ... ,b n ]- 

za~ lr P[l, b 2 - 1 ,... ,6„] 



+zV[l,b 2 , ...,b n ] 

, for b 2 > 2. 

V[l ,62,. 

• • ? b n ] 

= V[l, b 2 — 2,..., b n 

]-zV[l,b 2 -l,...,b n } 



+za b 2 V[b 3 ..., b n ) 

, for b 2 > 3 . 

V[bi,l,b 3 ,- 

• • ? b n ] 

= ^[6! -2,1,63,... 

, b„] - za}~ bl V[b 3 + 1,..., &„] 



+zV[ 61 - 1,1,63, 

• • •, b n ] for b 3 > 3 . 

P[ 2 ,l,& 3 ,. 

• ■ ? b n ] 

= aV[b 3 ,... ,b n ] - za 1 P[6 3 + 1,..., b n ] + zV[l, 1, 

V[l, 2,63,. 

■ • ? b n ] 

— P[b 3 + 1, • • •, b n ] - 

- zV[ 1,1, 63,..., b n ] + za 2 P[b 3 , 

^[l, 1 ,& 3 ,& 4 , • 

■ • ? b n ] 

= a~ b 3 V[l, 64 — 1,.. 

■, b n \ -zV[b 3 - t-l,& 4 ,...,6„] 



+zaV[b 3 l b i} ... ,6 

n] for 64 > 2. 

V[l, 1,63,1,65,. 

■ • ? b n ] 

= a b 3 V[l + 65,.. .,b n ] - zV[b 3 + 1,1,6 5 ,... ,b n ] 



+zaV[b 3 ,l,b 5 ,... 

? bn] • 


Proof. The statement follows from Lemma [ 3 ] and the discussion following it. □ 

Note that the standard braid-form diagrams appearing in both sides of any of the 
relations in Theorem [l] have the same parity. From the patterns of these relations, 
the second named author wrote a program in Mathematica® (see Appendix) which 
computes the Dubrovnik polynomial of any rational knot from a standard braid- 
form diagram. 

It is worth noting that although we began the reduction algorithm at the leftmost 
section of twists for programming purposes, beginning the reduction at the right 
hand side of the braid-form diagram reduces the number of cases needed to be 
considered. This is because the Dubrovnik skein relation, when applied to the right 
hand side of the diagram, does not result in diagrams (states) which are not in 
braid form, and thus does not require the sliding move. 


4. Coefficient polynomials and a reduction formula 

In this section we use the mechanics of the Dubrovnik skein relation described 
in Section [ 3 j to create an expression for the Dubrovnik polynomial of a standard 
braid-form diagram relative to the number of half-twists in a particular section. 

Our approach is motivated by the consistent recurrence of what we will call the 
coefficient polynomials of the Dubrovnik polynomial for the reduced diagrams. 

To begin with, we borrow a notation from [ 3 ] and consider a family of links L m 
(where m is an integer) which are identical except within a certain ball, where they 
have the segment indicated as in Figure [ 111 ] Thus we are considering links that are 
identical except for a chosen section of twists (see Figure 14 for an example). 
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L m , TO > 0 



L 0 L_ 



Lmi HI 0 



DC 


L+ 


L 


OO 


Figure 13 . The family of links L m 


For our purposes, L_ (or L + ) is used in conjunction with L m for to > 0 (or 
to < 0 ). We will explicitly show the case of D\b\,... , 6„] for all bi s positive (the 
negative case is treated similarly). The diagram L + or L_ appears in our reduction 
formula because of the possibility that we could have a section with one half¬ 
twist (unlike in [ 3 j, where all bi s are even). As discussed previously, applying the 
Dubrovnik skein relation to a section of one half-twist results in changes to the next 
section of twists, which L + and T_ address. 

We begin by considering a section of to > 0 half-twists in the rightmost block of 
twists in a standard diagram D[bi, ..., b n \, where n is odd and all b^s are positive. 
Figure fl 4 ] shows a family L m for to = 3 . 



Figure 14 . An example for L m , m = 3 


Remark 1 . In Section [ 3 ] we discussed how applying the Dubrovnik skein rela¬ 
tion affects a diagram. Observe that since the A-state reduces a diagram by one 
half-twist, it is equivalent to L m _ i, i.e., it has to — 1 half-twists in the chosen 
section. In addition, the switched-crossing state reduces the diagram by two half¬ 
twists, and therefore, it is equivalent to L m _ 2 . The trickier part is the relationship 
of the B-state with L m . The Dubrovnik skein relation applied to the rightmost 
crossing in a standard diagram D[bi,...,b n ] results in a B-state whose polyno¬ 
mial evaluation is a 1 ~ bn V[bi, ..., b n _ 1] for n odd and a bn ~ 1 V\bi 1 ..., & n _i] for n 
even. On the other hand, the polynomial evaluation, P(L 00 ), of the diagram is 
a m ~ brl 'P[bi, ..., 1] if n is odd and a bn ~ m V[bi 1 ..., b n - 1] if n is even. Therefore 

the evaluation of the B-state is a 1 ~ m P(L 00 ) for n odd and a m ~ 1 P(L 00 ) for n even. 
Therefore, the Dubrovnik skein relation applied at the rightmost crossing in the 
diagram D[b \,..., b n \ with n odd can be rewritten as 

( 4 . 1 ) P(L m ) = P{L m _ 2 ) - za 1 ~ m P(L 00 ) + zP{L m _ 1). 
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4.1. Recurrence relations for the coefficient polynomials. Let A rn denote 
the coefficient of P(L_), B rn the coefficient of P(L 0 ), and C m the coefficient of 
P(L oc ). We will show that we have the following recurrence relations: 


A ra — zA m _i A rn _ 2, Pm zPm—1 f Pm—2) and Cm — CL Cm—l 4” zPm—li 


where m > 2. Before that, we prove three lemmas involving these coefficients. 


LfJ 


Lemma 4. Let P m {z) = z m 2l \ ), for m £ N. Then 


i =o 


Pm — zPm — 1 4~ Pm — 2 i foi all 771 ^ 2. 


Proof. Let P m be as above, where m > 2. Then, we have: 


zPm—l f Pm—2 — Z 'y ' 


Z 2^ 2 

*=0 

L^J 

= E * 

2=0 

L^J 

= E * 


— 1 —2i/ TO ~ !-* 


1 -\ L?J 

m—2i I m ~ 1 - M + ^ 
i= 1 


E 

i=0 


v m—2i 


m — 2 — i 


1 A L?J 
ro-2i/ m_1 ~M 

i=l 


m — 1 — i 
i — 1 


m — 1 — i 
i — 1 


■ 2 . 


Case 1 : to odd. Then = LtJ ■ Thus we have 


L 2 J z m- 2 i z m- 2 i (m-l-i 

i— 1 


E 

2=1 

L¥J 

E* 

*=i 

LfJ 

y r m-2i (m-l 
i=1 

LfJ 

„m-2i fm.-I 


2=1 

ra — 1 — A (m — l — i 

L M i-i 




2=0 


Brn.- 
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Case 2 : to even. Then Lfy^J = LfJ — 1 > and we have: 


L^J 

E * 

i =1 

L¥J-i 


m—2i 


m — 1 — * 


LfJ 

E 

i=1 


y m—2i 


m — 1 — i 
i — 1 


= E 


^m—2i 


i—1 


m — 1 — i\ ( to — 1 — i 

i M <-i 


Pz n 

LfJ 

= E 


y m- 2L?j ( m - LfJ - 1 

1 LfJ -1 

L 2 J y m—2i (™ - A 2L® j f TO - LfJ - 1 

' l ■ i ‘I LfJ -1 

LttJ" 1 / \ 

= ^2 z m - 2i ( m M +-Z ™-2l?J 

i=0 \ 1 J 


m- LfJ - 1 
LfJ - 1 


Li 


= E 


i=0 


m—Zi I m 1 1 n 

= Dm. 


Therefore, the statement holds for both to odd and even. 


□ 


Lemma 5. Let A m (z) = ^ z m 1 2i ^ , for m € N. Then 

A m = T j4 m _2, for to 2 . 

Proof. Note that Ai = 1 and A rn = B m _ i for m > 1 , where B m is defined as 
in Lemma [ 4 ] Therefore, the statement follows by substituting to — 1 for to in 
Lemma [ 4 ] □ 


m L 1 ! 1 -! 

Lemma 6. Let C m (z, a) = ^ z-’~ 2 l a-’~ m 

2=1 *=o 

for to € N. Then the following hold: 



LfJ 

and B m (z) = ^ z m-2i 
»=o 



— a Cm— l T zB m — i, for to > 1 
c m = Cm—2 + za 1 ~ m + zC m - 1, for TO > 2 . 


Proof. Considering C m and as above, we have the following: 

m-lL^J 

= a” 1 ^ £ ^- 2 V- m+1 

2 = 1 j=0 

L^J 


i -1 - i 

i 


+* E 

2—0 

m L^2^J 


„m— 1—2i 


m — 1 — i 
i 


j-2i n j-m 


= EE 2 ' "“ a ' 
2=1 »=0 


-m /O' - 1 - * 


= C m . 
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Therefore, the first equality holds. We verify now the second equality. We observe 
first that 


Cm -2 + za 1 m + zC m -1 


2-2^1 


E E 


z j-2i a j-m+2 



+ za 1 


m 


m—l L^J 

+ Z e E ^“ 2 v- m+i 

1 = 1 2—0 

In addition, by dividing out by a~ m , we obtain 



Cm 

a -m 


m LS-J 

+£ £ 

1=2 i=0 



and 


C m -2 + za 1 T zCm—i 

a~ m 


m-2L i r t J 

E E 


j=i i=0 


—2i j+2 

z- 7 a J 



+ za 


m—l L^J 

+2 ^ ^ 0- 7_2 *a' 7+1 
j=i i=o 



Thus it suffices to show that 


P- 2 > £ £ 

1=2 2=0 


z j ~ 2 i a j 


j-l-i 
i 


m—l L 2 


= E E 

j = 1 2=0 


^- 2 V+ 2 


j - 1 - * 
i 


m— 1 L 2 


+ 2 E e 

1 = 1 2=0 


^- 2 V +1 


j ~ 1 - i 


Let’s take a look at the right hand side of the above equality. 


j- 2 i f.i - 1 - ' 


m —2 LtJ 

= £» )+2 £2. 

3 =1 *=0 

I J ~ 3 I 

m L 2 J / • _ q 

= J2 aj Y, zj ~ 2 ~ 2i J 1 

j =3 

I J ~ 3 I 
m L 2 -I 

= E « 3 E * 

1=3 2=0 

"L 1 ? 5 ! 


m—l LW£ 

E « J+1 E - 

j =1 i =0 

I J ~ 2 I 
m L 2 J 


j—21 f.i- 1 ' 


+ 2 


E« J E 


z- 7 


— l— 2 if J ^ * 


1=2 2=0 

m L^J 

xy e ^ 

1=3 2=0 


2-2-22 /O' - 3 - *\ ST J y J .2-22 f J - 2 - ^ + a 2 0 2 


= a-' 0 

1=3 


E 

*=o 


— 2 — 21 /0 ^ * 


L£mJ 

E 

i=0 


-2 i (J ^ ^ 


, 2 2 
H - a Z . 





A FORMULA FOR THE DUBROVNIK POLYNOMIAL OF RATIONAL KNOTS 


13 


Let S be the inside of the sum above. Then, we have 


5 = L £ \-2-2id-3-i 

i=0 

L^J 


E 

2=0 


-2 i lJ ^ i 


= E 




-2 i(j ~ 2 ~ i 


i — 1 


L^J 

E 


— 2 * (3 2 * 


i—1 ' ' i—1 

Case 1: If j is even then = LE^J- Therefore, we have 


+ 1 . 


L^J 

« = £ 

i— 1 

L^j 


i - 2 - i\ , (j-2-i 
i—1 


+ 1 


—2i l 3 1 ^ 


+ 1 


= E * 

i=1 

L^J 

= E 

2=0 


Case 2: If j is odd then \_^-\ = |_|j • Then, we have 


—2ifj 1 * 


L^J 


^ = E 


—2?: /0 2 i 

i — 1 


Lij-i 

E 




—2i(j 2 * 


i=l N 7 i=l 

Breaking off the last term of the first part, we get 


+ 1 . 


Rfcl-i 


^ = E 


-2i a - 2 - i 

i—1 


LiJ-i 

E 


-2*/j' 2 i 


+ 1 


+ 2 ;' 


-2[i^l(3 2 L J 2 1 J 

L^J -1 


LVJ-i 


- E 


-2 i 


2=1 

- 2 L^J 


, (j-2-i 

i—1 


+ 1 


+2 
L^J-i 


i-2-L^J 


: E 


— 2ifj 1 ^ 


2=1 


,- 2 L^J 


i- 2 -L^ 

L^J-i 


Since j is odd, j = 2p + 1 for some p € N. Then = Pi 


so 


j-2-L^J 
L^J-i 


j -p-2\ _ fj-p-1 
P-1 


= 1 , 


for the choice of j above. 
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Therefore, in either case, we get that the sum S satisfies the following equality: 




s= £ 


—2i (j 


i =0 


Hence, our full original equation becomes: 

Wj 


RHS = 


£< 

0 =3 


£ 


- 2-2 i (j ^ i 


i =0 

m L^J 

£« J ^ £ 2 


LtwJ 

£ 

i =0 


-2ifj ~ 2 - i 


, 2 2 
+ a z 


-2» / i - 1 - * 


4=3 


t=0 

L^J 


+ a 2 2 2 


^ z" 2l ' J 1 * 


4=2 


i=0 


which shows that the identity (4.2 1 holds. Consequently, the second equality in the 
statement holds. □ 

4.2. The reduction formulas. Lemmas [4] [ 6 ] and our previous analysis lead to 
the following result, which we call the first reduction formula (compare with jlSI , 
Theorem 1.14]; we will give more details on this comparison in Remark [2]). 

Theorem 2. Let m be a positive integer and consider a family of link diagrams 
L 0 , L_, Loo and L m which are identical except within a ball where they differ as 
indicated in Figure [73] Then 


(4.3) 

where 


P(L m ) = A m P(L_) + B m P(L 0 )-C m P(L oo ), 


■'177 

B n 

Crr 


= B m ^i for m > 1 and Ai = 1 

LfJ 

_ ^ z m ~ %i ( m ~ i 


i -0 


L^J 


= £ £ “" o; 




-m (j ~ 1 - * 


4 = 1 i= 0 

Proof. We will proceed by induction on to. If m = 1, then A\ = 1 , B\ = C± = z 
and the equation (4.31 is merely the Dubrovnik skein relation. Thus the statement 
holds trivially for to = 1 . 

Let to = 2. Then A 2 = z, B 2 = z 2 + 1, and C 2 = za~ l + z 2 . By the Dubrovnik 
skein relation, P(P 2 ) = P(L 0 ) — za~ 1 P(L 00 ) + zP(Li) and P{L\) = P(P_) — 
zP(Loo) + zP{Lq). Substituting -P(Li) in the identity for P(P 2 ), we obtain 

P(L 2 ) = P(Lo) — za~ 1 P(L 00 ) + zP(L_) — 0 2 P(P oo ) + z 2 P(L 0 ) 

= zP(L-) + ( z 2 + l)P(Lo) — i za 1 H" z 2 )P(Loo). 

Therefore, the equation (I4.3I) holds for to = 2. 
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Suppose that (4.3) holds for all 3 < i < m — 1, i £ N. In Remark [I] we showed 
that the Dubrovnik skein relation applied at the rightmost crossing in a standard 
braid-form diagram of odd length can be rewritten as follows: 

P(L m ) = P(i m _ 2 ) - za 1 ~ m P(L 00 ) + zPiLn-i). 


Note that by the way we defined L m , the diagrams Lq,L_ and Loo are the same 
regardless of the chosen m. Substituting our assumed formula in for P(L m _i) and 
P(i m _ 2 ), yields 

P{L m ) = [T m _ 2 P(L_) + P m _ 2 P(L 0 ) - C' m _ 2 P(P 00 )] - 2 a 1 - TO P(i 00 ) 

Az[A m -lP(L-) + P m _iP(Lo) — Cm-lP{Laoj\ 

= [Am —2 + zA m _i\P(L_) + [P m _ 2 + zB m _i}P(L 0 ) 

— [C m - 2 + za 1 m + zCm-llPi.Lao). 

Employing Lemmas [4] [6] we obtain 

P(L m ) = A m P(L_) + B m P(L 0 ) - C m P(L 00 ), 


which completes the proof. 


□ 


Remark 2. The coefficients A m {z) and B m (z) used in Theorem [2] are a version 
of the Chebyshev polynomials. (We thank Sergei Chmutov and Jozef Przytycki for 
pointing this to us.) The Chebyshev polynomials of the second kind, U m (x), satisfy 
the initial conditions 

Uq(x) = 1, U\(x) = 2x 

and the recursive relation 


U m {x) = 2xU m -l{x) - U m - 2 (x), for 771 >1. 

It is an easy exercise to verify that 

A m (z) = i 1 ~ m U m - 1 > f° r 171 — 1’ where i 2 = — 1 . 

We remark that Jozef Przytycki obtained (see m) a similar reduction formula as 
the one we gave in Theorem [2] (we thank him for telling us about this). The for¬ 
mula in Jl5l Theorem 1.14] presents the 2-variable Kauffman polynomial (not the 
Dubrovnik polynomial) of a link diagram L m with m > 0 in terms of the polyno¬ 
mials of the associated link diagrams L 0 ;^oo and L + (not L_). The polynomials 
v[ m \z) in m Theorem 1.14] are related to our coefficient polynomials A m (z ) as 
follows: 

A m (z) = i 1 ~ m v[ m \iz) for m > 1 . 


We consider now the case m < 0, to obtain the second reduction formula. Given a 
standard diagram D[bi,... ,b n ] with n even and &j > 0 for all 1 < * < n, we consider 
a block of to half-twists in the rightmost section of twists in the diagram. Since 
n is even and all bfs are positive, the rightmost section of twists in the diagram 
is in the upper row, which corresponds to the case to < 0. We give the resulting 
statement below. 


Theorem 3. Let m be a negative integer and consider a family of link diagrams 
To, L + , Loo and L m which are identical except within a ball where they differ as 
indicated in Figure [73] Then the following equality holds: 

(4.4) P(L m ) = A m P(L + ) + B m P(L 0 ) — C m P(P 00 ), 
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where 


A rn — 
B m = 

C m = 


B m +1 for m < —1 and A_i = 1 
L Djij 

{-z)^- 2x i 


2=0 


m\ — i 
i 


M L^J / . - . 

H 3 1 

j —1 i=0 


Proof. The proof is similar to the proof for Theorem [2j we use induction on \m\. If 
to = —1 then A-i = 1, P_i = C_i = — z and if to = —2 then A_ 2 = — z, P_ 2 = 
z 2 + 1, C- 2 = — za + z 2 . 

By the Dubrovnik skein relation, P{L_ 1 ) = P(L + ) — zP(L 0 ) + xP(L 00 ), and 
thus the formula (4.4) holds for to = —1. Moreover, P(i_ 2 ) = P(i 0 ) + zaP(Loo) — 
zP(L_ 1 ) and equivalently, 


P(P_ 2 ) = -zP(L+) + (z 2 + l)P(Lo) - (-za + z 2 )P(P 00 ). 


Thus the statement holds for to = —2. Now we suppose that the statement holds for 
all negative integers larger than to, and we prove it holds for to. By the Dubrovnik 
skein relation, we have 

P(L m ) = P(L m+2 ) + za^ m ~ 1 P(L 00 ) - zP(L m+1 ). 

A word of clarification is needed here, since to is negative. In the above notation, 
the diagrams L m + 1 and L m +2 contain 1 and, respectively, 2 less half-twists is the 
considered section of twists. By the induction hypothesis, 

P{L m + 1 ) = A m+ 1 P(L + ) + P m+ 1 P(L 0 ) — Cm+i-f > (ioo) 

P(I / m+ 2 ) = A m+ 2 P(L_|_) + B m+ 2 P(Lo) — Cm+ 2 P(Too)- 

Substituting these into the above expression for P(L m ), we obtain: 

P(L m ) = (A m+ 2 — zA m+ i)P(L + ) + (P m+2 — 2 iP m+ i)P(Lo) 

-(Cm +2 - za -™- 1 - zC' m+ 1 )P(L 00 ). 

Similar proofs as in Lemmas EMU can be given to show that for to < — 2 , we 
have 


A m = zA m+ i + A m+2 

Bm = zB m + 1 T B m + 2 

Cm. — zC^Tj-f-l 


— m +1 




m+2 > 


and therefore, the desired formula (4.41 holds for to. 


□ 


Remark 3. Consider a standard braid-form diagram D[bi ,..., b n ] with all 6 ,’s 
positive (and with n odd or even). Applying Theorems [2] and [3] recursively for to = 
b n , we can write the Dubrovnik polynomial P[b\,... ,b n ] in terms of polynomials 
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associated to standard braid-form diagrams with fewer sections of twists, as follows: 

L%-J 

V[b 1 ,...,b n ] = Y,(e n z ) b -- 2i ( bn 7 l )v[h,..., b n - 1 ,0] 

i =0 ' ' 

[NyAj 

i =0 ' ' 

b n L^J 

-EE (^r 2 v-f-‘"' 

7=1 i=0 


j ~ 1 - * 
i 


V[bi ,... oc], 


where e n = (—l) n U Using that U[ 6 i,..., 6 „_i, 0] = a £ " b " _ 1 P[ 6 i,..., 6 n _ 2 ] and 
that the following diagrams are equivalent as links: 


D[bi,.. -1] = £>[6 i,.. .,b n -1 - 1], D[bi,. .. ,&„_i,oo] = D[bi ,... 

we obtain that 

V[bi, ...,b n ] = a '" 6 - 1 ( K 7 *W i,. • •, b n _ 2 ] 

i =0 \ 1 / 


L^J 


(4.5) 


+ X] 1 

i=0 

6 „ L^J 

-X! X! -” a 

7=1 »=0 


Vi 6 !.'”. 6 "-!- 1 ] 


j-2i n e n (j-b n ) (J 1 * 


U[6i,..., 6 n _i], 


for all n > 2. Note that the sign e n = (—l ) n_1 is necessary for switching in 
our computations between m < 0 and m > 0 (see Figure 131, as the calculations 
switch between reducing lower and upper rows of a standard braid-form diagram. 
Observe that reducing a lower row (or upper row) corresponds to an odd-length (or 
even-length) standard braid-form diagram. 

For a given standard diagram D[b\, ..., b n -\, b n \ with all bj/s positive and n > 2 
(where n is odd or even), let 

Xn,b n ■= V[bi, ■ ■ • i,b n \ and a; ni6 „_i := V\bi _,... , 6 n _i, 6 n - 1 ]. 


More generally, let 


x n ,k :=V[bi,... ,b n -!,k], where k < b n 


and 

®n-i, 6 „-i := P[&i,..., b n - 1 ] and x n - 2 ,b n -2 ■= P[bi, ..., 2 ]- 

Thus the first subscript n in x n ^ stands for the length of a standard diagram and 
the second subscript k corresponds to the number of half-twists in the rightmost 
section of twists in the given diagram. It is important to note that the integers 
61, ..., b n _ 1 are fixed in the notation for x n ^. 
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With this new notation, the formula in (4.5) becomes: 


l¥j 


x n ,b n = a ( 1)71 b " _1 E ((“!)" ±z ) 


-l~\b n -2ilbn * 


i=0 


%n—2,b n —2 


L^J 


(4.6) 


i-l _\i>„-l-2i I bn 1 * 


%n—l,b n —i — 1 

i - 1 - i 


+ E ((- 1 )” -1 ^) 6 

2=0 

-E E ((- 1 )"- 1 ^'- 2< a ( - 1)B “ 1 0-- 6 » ) r , 

j=i *=o 

This formula gives the Dubrovnik polynomial of a standard braid-form diagram of 
length n in terms of the polynomials of standard diagrams of lengths n — 1 and 
n — 2, greatly reducing the amount of the necessary iterations, when compared to 
the defining skein relation alone. 

For a standard diagram D[bi,, 6„_i, b„\ with 6,'s positive and n > 2, we let 


(4-7) r n,b„ 


= a^ n ~ lb - 


L 2 J 

. £((-i 

2 = 0 ' ' 


i-l ~\b n -l-2i I t>n 1 * 


[ b n-l j 

(4-8) P„, 6n : = E (M)" -1 *) 

2=0 

6 n L^J 

(4.9) Z n , 6 „: = -EE((- i r lz ) i " 2iflHrlH) 

i=i i=0 


j ~ 1 - * 
i 


The identity (4.6) becomes 


Xn,b n ^n^bn^n—l.bn — i 4 Pn,6 n *^n-l,6„_i-l 

or equivalently, 

(4.10) 

— ln,b n P [6i, . . . , 6 n _i] C 7 ^ [6‘| , . . . , 6^—2] T Pn,b n 7 * [61, ... , b n — 1 1 ]. 


Note that the terms in equation (4.10) depend on the fixed integers 6 ra _i and 


b n — 2 : which are determined by the form of the original standard diagram D[b\,... ,b n . 

Applying the Dubrovnik skein relation for diagrams D[bi] and £>[ 61 , 62 ] with 
both 61 and 62 positive, we see that the identities (4.6) and (4.10) hold for n = 1 
and n = 2 if we set 

xo,b 0 = 1, £-i,b_i = z~ 1 a + 1 - z~ 1 a~ 1 and x 0 ,b o -i = a -1 - 
We define 60 = 0 and 6 _i =0, and therefore 

x 0 ,o 4 

X—ib_ 1 = a.’-i,o = z~ 1 a + 1 — z~ 1 a~ 1 
Xo,bo — l *»o,—1 a 

= 1 . 


Note that xi,o = P ((E) 


Therefore, with the above conventions, the equalities (4.6) and (4.10) hold for 
all n £ N. 
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5. A CLOSED-FORM FORMULA FOR THE DUBROVNIK POLYNOMIAL 


Now that we have a formula for the Dubrovnik polynomial of a standard braid- 
form diagram in terms of polynomials of standard diagrams of shorter lengths, it 
is one more step to obtain a closed-form expression for the Dubrovnik polynomial 
of the original diagram. Specifically, we seek a closed-form expression for x n: b n 
in terms of a, z and the entries of the n-tuple (&i, ..., 6„_i, b n ) associated with a 
standard diagram D[b\,, b n _ i, b n \. For this, it is convenient to consider first a 
more general situation represented by the following lemma. 


Lemma 7. Given a standard braid-form diagram D[bi ,..., b n -\, b n \ of length n 
and with bi > 0 for all 1 < * < n, let Lj, rij and pij be elements of a certain com¬ 
mutative ring R (here we are interested in R = Z[a ±:L , z ±:L ]) and define recurrently 
the sequence x nt b„,n > —1, of elements in R by the relation 

(5.1) Xn,b n — ln,b n 'Xn—l,b rL —\ T X n ^b n Xn—2,b n —2 T Pn,b n 'En—l,b Tl —i — li fbl Tl P 1 


where iC-po, ^o.cu and £ 0,-1 are fixed elements in R , and bo = 5_i = 0. Let F 
be the set of all strictly decreasing integer sequences / = {/i, fi, ■ ■ ■ , /;}, where 
fi = n, fi — fi + 1 = 1 or 2, fi = 0 or — 1, and only one of 0 or —1 is present in a 
sequence / (that is, if fi = —1 then /;_i ^ 0). Each fi = Ci or di , and if fi = di 


then fi - f i+ 1 = 1. For each f t £ /, let t fi = < 

[0 if fi _i = c.i_i. 

Then x nt b n can be written in terms of the initial conditions ai-po, Xo,o, and £o,-i 
and the elements hj, r^j and Pij as follows: 


(5.2) x n ^ n 


where 


Ifi,bf i —tf i r f i ,b fi -t fi Pf i ,bf i —t/ i 
*6A(/) i6p(/) *S7(/) 


Hf) = {i\fi = c i and a - f i+1 = 1 } 

p{f) = {*| fi = Ci and d - fi+ 1 = 2} 

7 (/) = {i\fi = ^ and di - f i+1 = 1 }. 


Proof. The formula describes the tree of calculations resulting from successive ap¬ 
plications of the recurrence formula We represent the computational tree 

using a layered rooted tree (as in w, where each layer i collects the coefficients 
li.ji r i.j i and pij. Moreover, each layer i corresponds to all terms Xij whose first 
subscript is i. The Ledges and r-edges are drawn on the left and middle, respec¬ 
tively, at a vertex in the tree. According to the formula for x n ^ n , the Ledges are 
of length 1 and r-edges of length 2; that is, the Ledges connect vertices in the tree 
that are at one level apart, and the r-edges connect vertices that are at two levels 
apart. Similarly, the p-edges are of length 1. We draw the p-edges to the right, 
and color them red. Figure |T5| shows a computational tree with n = 3. The set of 
integer sequences F generates the possible paths from the root to a leaf of the tree. 
Each sequence / corresponds to a unique branch in the tree, with fi representing a 
vertex of the path. An entry Ci in a sequence / corresponds to traveling along the 
left or middle edge incident to a vertex /, located at the Lth level in the tree. On 
the other hand, an entry di in / corresponds to traveling along the right edge at 
that vertex. Any path in the tree produces a product of L, r-, and p-coefEcients, 
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3 

2 

1 

0 

-1 



Figure 15. Computational tree for n = 3 


and the expression for x n ^„ is the sum over all paths in the tree, where each path 
starts from the highest level and ends at a leaf of the tree. □ 

To obtain our desired closed-form formula for the Dubrovnik polynomial of a 
standard braid-form diagram D[bi,... ,b n \ with all bi > 0 positive, we combine 
Lemma [7] with the reduction formulas provided in Theorem [ 2 ] and Theorem [ 3 ] (and 
the associated formula ( 4.10[ ) in Remark [ 3 ]). Specifically, we associate the vertices 
in the layered tree of Lemma [T] with the (Dubrovnik polynomial of) rational knot 
diagrams in standard braid-form obtained by successive applications of the reduc¬ 
tion formulas in Theorem [2] and Theorem [3j The root of the tree corresponds to the 
given rational knot diagram. Each vertex v in the tree is associated with a diagram 
L m as in Theorem[2] (if m > 0) or Theorem[3] (if m < 0), and the three descendants 
of v correspond to the diagrams L^, Lq and L_ (or L + ) associated with L m . Figure 
[IS] depicts the computational tree for a standard braid-form diagram of length 3. 

Notice that the /-, r-, and p-coefficients correspond to the coefficient polynomials 
C m (z,a), B m (z) and A m (z), respectively, introduced in Section |4j The p-edges in 
a computational tree are highlighted in red to differentiate them from the other 
edges in the tree, as they reduce a different standard diagram than the other edges, 
namely D[b \,..., bi — 1], decreasing not only the number of sections of twists in the 
diagram but also the number of half-twists in the rightmost section of twists. 


Theorem 4. 

all 1 < i < n. 


Let D[bi, 
Then 


, b n \ be a standard braid-form diagram with bi > 0 for 


V[bi = 


bfi-tfi Lbrd 


n 

E E 

/£F iex(f) 

j— 1 *=0 



■ n al - 1 )* -1 **- 1 


i6p(/) 

2=0 

b f. —t f. — 1 

L J 

n e ((- 

-1 


~ 1 zV~ 2i a ( -~ 1 '> U ti~ b fi +t fi) 


j -1-i 


where 


i=0 

• F is the set of all strictly decreasing integer sequences f = {/ 1 , / 2 , • • •, fi}, 
where fi = n, fi — fi + 1 = 1 or 2, fi = 0 or —1, and only one of 0 or —1 is 
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present in the sequence f (that is, if fi = — 1 then fi-\ ^ 0). Each fi = Ci or 

di, and if fi = di then fi-fi+i = 1. For fi £ f,t fi = •P ^ ^ _1 

[0 if fi -1 = Ci- 1 . 

• Hf) = {*1 fi = Ci and Ci - fi+i = 1 }, 

• P(/) = {*1 fi = Ci and Ci - f i+1 = 2}, 

• 7 (/) = {i\fi = di and di - f i+1 = 1 }, 

• *o,o = 1 ) £-i,o = z~ x a + 1 — z~ 1 a~ 1 and Xo,-i = a 1 - 

Proof. The statement follows from Lemma [7] and Remark [3j □ 


The Dubrovnik polynomial of a standard braid-form diagram D[bi,, b n ] with 
all b t < 0 is obtained from the closed-form formula (given in Theorem [4]) for the 
Dubrovnik polynomial for the diagram D[|&i|£> n |], in which one is applying 
the replacements a <—> a^ 1 and z < — > —z. That is, 

V[bi ,..., b n ](a, z) = V[\h\,..., |&„|](a -1 , -z), h <0, VI < i < n. 


Example 1. We compute the Dubrovnik polynomial for the standard diagram 


D[A, 3, 5] depicted in Figure 16 



Figure 16. Standard diagram D{ 4,3,5] 


The length of the diagram is n = 3 and there are 17 possibilities for the sequences 
in F, which are given in the first column of Table |T] To differentiate cfs and 
df s in a sequences / £ F, we denote df s with a subscript p. In this table, the 
second and third columns list the terms c» £ / corresponding to i £ A (/) and 
i £ p(f), respectively. The fourth column lists the terms di £ / corresponding to 
i £ 'y(f). The last column collects the product of the coefficients corresponding to 
each branch in the tree; each such product is a term in the Dubrovnik polynomial 
V[A, 3, 5] associated with the standard diagram D[ 4,3, 5]. 

The tree that corresponds to a standard diagram of length 3 is given in Figure [T5| 
Recall that each sequence / corresponds to a path in the tree starting from the 
root at level 3 and ending at a leaf, recording the level for each vertex in the path 
and writing a subscript p (which marks df s) for the vertices where we chose the 
rightmost edge highlighted in red. 

To see how and why the computational tree works, let the vertices of the tree 
correspond to the following polynomials, which we evaluate using the formula (5.1) 
in Lemma [7j 


ms, 5] 

m 3] 
m 2] 
m 
m 


i3,5P[4,3]+r 3 ,5P[4]+P3, s P[4,2] 
+ T 2,3^0,0 +P2,37 > [3] 

h,2P\A\ + T 2 , 2 £ 0,0 + P2,fP[ 3] 
^ 1 , 4 £ 0,0 + Tl, 4 £—1,0 + Pi, 4 ^ 0,-1 

^1,3*0,0 + r l,3*-l,0 +7 , 1,3*0,-1- 
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/ 

Hf) -+ {+} 

P{f) -+ {+} 

7(/) -+ {<*»} 

products of coefficients 

{3,2,1,0} 

{3,2,1} 

0 

0 

^ 3 , 5 ^ 2 , 3 ^ 1 , 4 * 0,0 

{3,2,1,-1} 

{3,2} 

{1} 

0 

^ 3 , 5 ^ 2 , 321 , 4 *— 1,0 

{3,2, l p , 0} 

{3,2} 

0 

{l P } 

h,bh,3Pl, 4 X 0-1 

{3,2,0} 

{3} 

{2} 

0 

^ 3 , 522 , 3 * 0,0 

{3,2 P , 1,0} 

{3,1} 

0 

{2p} 

h,5P2, 3h, 3*0,0 

{3,2 p , 1,-1} 

{3} 

{1} 

{2p} 

h,0P2,3 r 1,3*—1,0 

{3,2 P , l p , 0} 

{3} 

0 

{2p,lp} 

^3,5P2,3Pl,3*0,-l 

{3,1,0} 

{1} 

{3} 

0 

23,5^1,4*0,0 

{3,1,-1} 

0 

{3,1} 

0 

23,521,4* —1,0 

{3, l p , 0} 

0 

{3} 

{lp} 

23,54*1,4*0,-1 

{3 P , 2,1,0} 

{2,1} 

0 

{3p} 

P3,5^2,2^1,4*0,0 

{3 P , 2,1,-1} 

{2} 

{1} 

{3p} 

P3,5^2,2*1, 4 *—1,0 

{3 P , 2, l p , 0} 

{2} 

0 

{3p,Ip} 

P3,5^2,2Pl,4*0,-l 

{3p, 2,0} 

0 

{2} 

{3p} 

P3,5*2,2*0,0 

{3 P , 2 P , 1,0} 

{1} 

0 

{3p, 2p} 

4*3,54*2,2^1,3*0,0 

{3p, 2p, 1, —1} 

0 

{1} 

{3p, 2p} 

4*3,5P2,2*1, 3 * —1,0 

{3p,2p,lp,0} 

0 

0 

{3p,2p,lp} 

P3,5P2,2Pl,3*0,-l 


Table 1. Paths and terms for computing 5] 


Thus 

ms] = 


and 

= 


Therefore, we have 


ms, 5] = 


^2,3^14] + 22,3*0,0 + P2,3'P[3] 

h,3{h,4X0,0 + Tp42-1,0 + Pi,4*o,-i) + r 2 ,3X0,0 

+P 2 , 3(ll, 3X0,0 + Xl, 32 - 1,0 + Pi,3*0,-l) 

^2,3^1,4*0,0 + ^2,32l,4*-l,0 + ^2,34*1,4*0,-1 + r 2,3 x 0,0 
+P2,3h,3Xo,0 + P 2 , 32 l, 3 *- 1,0 + P2,3Pl,3Xo,-l 


^2,2'P[4] + 22,2*0,0 + P2,2^[3] 

h, 2 {li, 4 X 0,0 + 2i,4*-i,o + Pi, 4 X 0 - 1 ) + r 2 , 2 X 0,0 

+P 2 , 2 (^ 1 , 3 * 0,0 + 2i ; 3 2_i,o + Pi, 320 ,- 1 ) 

^ 2 , 2 ^ 1 , 4 * 0,0 + ^ 2 , 2 * 1 , 4 * 2 — 1,0 + l 2 , 2 Pl,4Xo,-l + 22 , 2 * 0,0 
+P 2 , 2h, 3*0,0 +P2,22l,3*-1,0 + P2,2Pl,3*0,-l* 


; 3 ,5^[4,3] + 2 3 , 5 P(*D[ 4 ]) +p 3j5 P[4,2] 

^3,5^2,3^1,4*0,0 + ^3,5^2,321,4*— 1,0 + ^3,5^2,3Pl,4*0,-l 
+^3,522,3*0,0 + ^3,5P2,3^1,3*0,0 + ^3,5P2,32l,3*-l,0 
+^3,54*2,34*1,3*0,-1 + *3,5^1,4*0,0 + 23 , 52442 - 1,0 
+2 3 ,5Pl,4*0,-l + P3,5^2, 2 ^ 1 ,4*0,0 + P3,5^2,22l,4*-1,0 
+P3,5^2,2Pl,4*0,-l +P3,5*2, 2 * 0,0 + P3.5P2, 2 ^ 1 ,3*0,0 
+P3,5P2, 2 * 1 , 3 2_i, 0 +P3,5P2,2Pl, 3*0, — 1 * 
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Note that this result coincides with the sum of the products in the last column 
in Table [1] which are obtained from the computational tree in Figure [T5| and its 
associated paths. Using the relations (4.71 through (4.9) and the chart of paths, 
we obtain 


P[4,3,5] = 


^i + - i -2a 

, 2,2 


o 4 20 Z _ , 3 2z 2 3 z 2 6z 2 

- 3a 4 —H— — 5 az + 2 a z — 18 z H-„—I-^—I- -r- 

a 5 a a 8 a 6 a 4 


— 6 a 2 0 2 + 13 a 4 z 2 + + 


3z 3 12 ^ 3 50 3 5z 3 


a' 


a° 


a° 


5z 4 

~n6~ 


17 z 4 


1 C te¬ 


ar a 
,3,5 


a 

- 6 a 2 z 4 + 16 a 4 0 4 + —+ 
a' 


^ + —-26a^3 + a 3 z 3 - A3z 4 + ^ 


30 6 


-33 az 3 - 6a z - 360 b + —r + 

a° 


12z 6 12z 6 

+ 


10z 5 

, 2,6 


12 z 5 


15z 5 

a 

4 z 7 


a* 


a* 


+ 2a z + 7a z + - ir 
a° 


9 z 7 10 z 1 7 O 7 O 4z 8 7 z 8 2 8 4 8 

H—H-1 8az 7 - 5 a 3 z 7 - 16z 8 H-j- H-=- + 4a 2 z 8 + a 4 z 8 

a 6 a a 4 a z 


~ + — - 6 a 0 9 - a 3 0 9 - 3z 10 + 
a 6 a 


2z 


10 


+ a 2 z 10 + -- az 

a 


,n 


where we used that 


zo,o 

h,5 

T3,5 

h ,3 

h, 4 

T 1,4 
U,3 

h,2 


1, = x x a 


1 — 0 4 a 1 


-l 


li 


x 0 ,-i — a 

0 a -2 + 0 4 a _1 + 20 2 a ~ 1 + 0 5 + 3 z 3 + z) 
2 i ,4 


— 1 i ,2,-3 , ,3-2 

— —1 yza z a z a 

= a 3 (z 5 + 4z 3 + 3z), p 3 , 5 = 1 + 3z 2 + z 
= -1 (-z - za 2 + z 2 a - z 3 ), r 2j3 = - 1(20 + z 3 )a~ 4 , p 2j3 = 1 + 0 2 

= — l( 0 a -3 + z 2 a~ 2 + 0 3 a _1 + za ~ 1 + 0 4 + 20 2 ) 

= 0 4 + 30 2 + l, pi t i = z 3 + 2z 

= -(0 + zaT 2 + z 2 aT l + 0 3 ), r lj3 = 20 + 0 3 , p 13 = 1 + 0 2 


= 0 a — 0 


T 2,2 = a (0 +1), P2,2 = ~Z. 


6. Appendix 

This appendix contains the Mathematica® code written by the second named 
author, which computes the Dubrovnik polynomial V[b\, b 2 , • • •, b n ] based on suc¬ 
cessive iterations of the Dubrovnik skein relation applied to the left hand side of a 
standard braid-form diagram of a rational knot, as shown in Section [3] 

Clearff, a, c, d, bl, z, b2, b3, b4, b5, btail, i] 
f [1, {111 = c; 
f[l, {2}] = d; 
f[l, lbl_>] := 

f[l, {bl}] = f[l, {bl - 2 }] - z (a~(1 - bl) - f[1, {bl - 1}]) 


f[2, {1, 1}] := d; 

f [2, {1, 2}] := f [2, {1, 2}] = a~ (-1) - z (d - a‘2); 
f[2, {2, 1}] := f[2, {2, 1}] = a - z a~(-2) + z (f[2, {1, 1}]); 
f[2, {1, b2_}] := 

f[2, {1, b2}] = f[2, {1, b2 - 2}] - z (f[2, {1, b2 - 1}] - a~b2); 
f [2, {2, b2_>] := f[2, {2, b2>] = 
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a~b2 - z (a~(-l) f[2, {1, b2 - 1}] - f[2, {1, b2>]); 

f [2, {bl_, 1>] := f [2, -Cbl, 1}] = f[l, {bl + 1}]; 

f [2, {bl_, b2_}] := f[2, {bl, b2}] = 

f[2, {bl - 2, b2>] - z (a~(l - bl) f[2, {1, b2 - 1}] 

- f[2, {bl - 1, b2}]) 

f[3, {1, 1, b3_>] := f[3, {1, 1, b3>] = 

a~(-b3) - z (f[l, {b3 + 1>] - a f[l, {b3>]); 
f [3, {1, 2, b3_>] := f[3, {1, 2, b3>] = 

f[l, {b3 + 1>] - z (f[3, {1, 1, b3>] - a~2 f[l, {b3>]); 
f [3, {2, 1, b3_>] := f[3, {2, 1, b3>] = 

a f [1, {b3>] - z (a'(-l) f[l, {b3 + 1>] - f[3, {1, 1, b3>]); 
f [3, {1, b2_, b3_>] := f[3, {1, b2, b3>] = 

f[3, {1, b2 - 2, b3>] - z (f[3, {1, b2 - 1, b3>] 

- a~b2 f[l, {b3}]); 

f[3, {2, b2_, b3_>] := f[3, {2, b2, b3>] = 

a~b2 f[l, {b3>] - z (a~(-l) f[3, {1, b2 - 1, b3>] 

- f[3, {1, b2, b3>]); 

f[3, {bl_, 1, b3_>] := f[3, {bl, 1, b3>] = 

f[3, {bl - 2, 1, b3>] - z (a~(l - bl) f[l, {b3 + 1}] 

- f[3, {bl - 1, 1, b3>]); 
f [3, {bl_, b2_, b3_>] := f[3, {bl, b2, b3>] = 

f[3, {bl - 2, b2, b3>] - z (a~(-bl + 1) f[3, {1, b2 - 1, b3>] 

- f[3, {bl - 1, b2, b3>]) 

f[4, {1, 1, b3_, 1>] := f[4, {1, 1, b3, 1>] = 

a~(-1 - b3) - z (f [2, {b3 + 1, 1>] - a f[2, {b3, 1}]); 

f[4, {1, 1, b3_, b4_>] := f[4, {1, 1, b3, b4>] = 
a~(-b3) f[2, {1, b4 - 1>] - z (f[2, {b3 + 1, b4>] 

- a f[2, {b3, b4>]); 

f[4, {2, 1, b3_, b4_>] := f[4, {2, 1, b3, b4}] = 

a f[2, {b3, b4>] - z (a~(-l) f[2, {b3 + 1, b4>] 

- f[4, {1, 1, b3, b4}]); 

f[4, {1, 2, b3_, b4_>] := f[4, {1, 2, b3, b4>] = 
f[2, {b3 + 1, b4}] - z (f [4, {1, 1, b3, b4>] 

- a~2 f[2, {b3, b4}]); 

f[4, {1, b2_, b3_, b4_>] := f[4, {1, b2, b3, b4>] = 

f[4, {1, b2 - 2, b3, b4>] - z (f[4, {1, b2 - 1, b3, b4}] 

- a~b2 f[2, {b3, b4>]); 

f[4, {2, b2_, b3_, b4_>] := f[4, {2, b2, b3, b4>] = 

a~b2 f[2, {b3, b4>] - z (a~(-l) f[4, {1, b2 - 1, b3, b4}] 

- f[4, {1, b2, b3, b4>]); 

f[4, {bl_, 1, b3_, b4_>] := f [4, {bl, 1, b3, b4>] = 

f[4, {bl - 2, 1, b3, b4>] - z (a~(l - bl) f[2, {b3 + 1, b4>] 

- f[4, {bl - 1, 1, b3, b4>]); 

f[4, {bl_, b2_, b3_, b4_}] := f[4, {bl, b2, b3, b4>] = 

f [4, {bl - 2, b2, b3, b4>] 

- z (a~(-bl + 1) f[4, {1, b2 - 1, b3, b4}] 
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- f[4, {bl - 1, b2, b3, b4>]) 

f[5, {1, 1, b3_, 1, b5_>] := f[5, {1, 1, b3, 1, b5>] = 
a"(-b3) f [1, {1 + b5>] - z (f[3, {b3 + 1, 1, b5>] 

- a f[3, {b3, 1, b5>]); 

f[5, {1, 1, b3_, b4_, b5_>] := f[5, {1, 1, b3, b4, b5}] = 
a~(-b3) f[3, {1, b4 - 1, b5>] - z (f[3, {b3 + 1, b4, b5}] 

- a f[3, {b3, b4, b5}]); 

f[5, {2, 1, b3_, b4_, b5_>] := f[5, {2, 1, b3, b4, b5>] = 

a f[3, {b3, b4, b5>] - z (a~(-l) f[3, {b3 + 1, b4, b5>] 

- f[5, {1, 1, b3, b4, b5>]); 

f[5, {1, 2, b3_, b4_, b5_>] := f[5, {1, 2, b3, b4, b5}] = 

f[3, {b3 + 1, b4, b5>] - z (f[5, {1, 1, b3, b4, b5}] 

- a~2 f[3, {b3, b4, b5>]); 

f[5, {1, b2_, b3_, b4_, b5_>] := f[5, {1, b2, b3, b4, b5>] = 

f[5, {1, b2 - 2, b3, b4, b5>] - z (f[5, {1, b2 - 1, b3, b4, b5>] 

- a~b2 f[3, {b3, b4, b5>]); 

f[5, {2, b2_, b3_, b4_, b5_}] := f[5, {2, b2, b3, b4, b5}] = 

a~b2 f[3, {b3, b4, b5}] - z (a~(-l) f[5, {1, b2 - 1, b3, b4, b5>] 

- f[5, {1, b2, b3, b4, b5>]); 

f [5, {bl_, 1, b3_, b4_, b5_>] := f[5, {bl, 1, b3, b4, b5>] = 

f[5, {bl - 2, 1, b3, b4, b5>] 

- z (a~(1 - bl) f[3, {b3 + 1, b4, b5>] 

- f[5, {bl - 1, 1, b3, b4, b5}]); 

f [5, {bl_, b2_, b3_, b4_, b5_>] := f[5, {bl, b2, b3, b4, b5>] = 

f[5, {bl - 2, b2, b3, b4, b5>] 

- z (a~(-bl + 1) f[5, {1, b2 - 1, b3, b4, b5>] 

- f[5, {bl - 1, b2, b3, b4, b5>]) 

f[i_, {1, 1, b3_, 1, b5_, btail_}] := 

f[i, {1, 1, b3, 1, b5, btail}] = a~(-b3) f[i - 4, {1 + b5, btail}] 

- z (f[i - 2, {b3 + 1, 1, b5, btail}] 

- a f[i - 2, {b3, 1, b5, btail}]); 

f[i_, {1, 1, b3_, b4_, btail_}] := f[i, {1, 1, b3, b4, btail}] = 

a~(-b3) f[i - 2, {1, b4 - 1, btail}] 

- z (f[i - 2, {b3 + 1, b4, btail}] 

- a f[i - 2, {b3, b4, btail}]); 

f[i_, {2, 1, b3_, btail_}] := f[i, {2, 1, b3, btail}] = 

a f[i - 2, {b3, btail}] - z (a~(-l) f[i - 2, {b3 + 1, btail}] 

- f[i, {1, 1, b3, btail}]); 

f[i_, {1, 2, b3_, btail_}] := f[i, {1, 2, b3, btail}] = 

f[i - 2, {b3 + 1, btail}] - z (f[i, {1, 1, b3, btail}] 

- a~2 f[i - 2, {b3, btail}]); 

f[i_, {1, b2_, btail_}] := f[i, {1, b2, btail}] = 

f[i, {1, b2 - 2, btail}] - z (f[i, {1, b2 - 1, btail}] 

- a~b2 f[i - 2, {btail}]); 

f[i_, {2, b2_, btail_}] := f[i, {2, b2, btail}] = 

a~b2 f[i - 2, {btail}] 



26 


CARMEN CAPRAU AND KATHERINE URABE 


- z (a~(-l) f[i, {1, b2 - 1, btail}] - f[i, {1, b2, btail}]); 
f[i_, {bl_, 1, b3_, btail }] := f[i, {bl, 1, b3, btail}] = 

f [i, {bl - 2, 1, b3, btail}] 

- z (a~(l - bl) f[i - 2, {b3 + 1, btail}] 

- f[i, {bl - 1, 1, b3, btail}]); 

f[i_, {bl_, b2_, btail_}] := f[i, {bl, b2, btail}] = 

f [i, {bl - 2, b2, btail}] 

- z (a~(-bl + 1) f[i, {1, b2 - 1, btail}] 

- f[i, {bl - 1, b2, btail}]) 

c : = a; 
e := a~(-l); 

d := z~(-l) a - z~(-l) a~(-l) + 1 - z a~(-l) + z a; 

The code for negative bi s is given below. 

Clear[f, a, c, d, bl, z, b2, b3, b4, b5, btail, i] 
f [1, {-1}] = e; 
f[l, {-2}] = d; 
f[l, {bl_}] := 

f[l, {bl}] = f[l, {bl + 2 }] + z (a~(-1 - bl) - f[1, {bl + 1}]) 
f[2, {-1, -1}] := d; 

f[2, {-1, -2}] := f [2, {-1, -2}] = a + z (d - a~(-2)); 
f[2, {-2, -1}] := f[2, {-2, -1}] = a~(-1) 

+ z (a~2 - f[2, {-1, -1}]); 

f[2, {-1, b2_}] := f[2, {-1, b2}] = 

f[2, {-1, b2 + 2}] + z (f [2, {-1, b2 + 1}] - a~b2); 

f[2, {-2, b2_}] := f[2, {-2, b2}] = 

a~b2 + z (a f[2, {-1, b2 + 1}] - f[2, {-1, b2}]); 

f[2, {bl_, -1}] := f[2, {bl, -1}] = f[1, {bl - 1}]; 

f[2, {bl_, b2_}] := f[2, {bl, b2}] = 

f[2, {bl + 2, b2}] + z (a“(-l - bl) f[2, {-1, b2 + 1}] 

- f[2, {bl + 1, b2}]) 

f[3, {-1, -1, b3_}] := f[3, {-1, -1, b3}] = 

a~(-b3) + z (f[l, {b3 - 1}] - a~(-l) f[l, {b3}]); 
f[3, {-1, -2, b3_}] := f[3, {-1, -2, b3}] = 

f[l, {b3 - 1}] + z (f [3, {-1, -1, b3}] - a~(-2) f[l, {b3}]); 
f[3, {-2, -1, b3_}] := f[3, {-2, -1, b3}] = 

a~(-l) f [1, {b3}] + z (a f[1, {b3 - 1}] - f[3, {-1, -1, b3}]); 

f[3, {-1, b2_, b3_}] := f[3, {-1, b2, b3}] = 

f [3, {-1, b2 + 2, b3}] + z (f[3, {-1, b2 + 1, b3}] 

- a~b2 f[l, {b3}]); 

f[3, {-2, b2_, b3_}] := f[3, {-2, b2, b3}] = 
a~b2 f [1, {b3}] + z (a f[3, {-1, b2 + 1, b3}] 

- f[3, {-1, b2, b3}]); 

f [3, {bl_, -1, b3_}] := f[3, {bl, -1, b3}] = 

f[3, {bl + 2, -1, b3}] + z (a~(-l - bl) f[1, {b3 - 1}] 

- f[3, {bl + 1, -1, b3}]); 

f [3, {bl_, b2_, b3_}] := f[3, {bl, b2, b3}] = 
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f[3, {bl + 2, b2, b3>] + z (a~(-bl - 1) f[3, {-1, b2 + 1, b3}] 

- f[3, {bl + 1, b2, b3>]) 

f[4, {-1, -1, b3_, -1>] := f[4, {-1, -1, b3, -1>] = 

a~(1 - b3) + z (f [2, {b3 - 1, -1>] - a'(-l) f[2, {b3, -1}]); 
f[4, {-1, -1, b3_, b4_>] := f[4, {-1, -1, b3, b4}] = 
a~(-b3) f[2, {-1, b4 + 1>] + z (f[2, {b3 - 1, b4>] 

- a~(-l) f[2, {b3, b4>]); 

f[4, {-2, -1, b3_, b4_>] := f[4, {-2, -1, b3, b4>] = 

a" (-1) f[2, {b3, b4>] + z (a f[2, {b3 - 1, b4>] 

- f[4, {-1, -1, b3, b4>]); 

f[4, {-1, -2, b3_, b4_>] := f[4, {-1, -2, b3, b4>] = 

f[2, {b3 - 1, b4>] + z (f[4, {-1, -1, b3, b4>] 

- a~(-2) f[2, {b3, b4>]); 

f[4, {-1, b2_, b3_, b4_>] := f[4, {-1, b2, b3, b4}] = 

f[4, {-1, b2 + 2, b3, b4>] 

+ z (f [4, {-1, b2 + 1, b3, b4>] - a~b2 f[2, {b3, b4>]); 
f [4, {-2, b2_, b3_, b4_>] := f[4, {-2, b2, b3, b4>] = 

a~b2 f[2, {b3, b4>] + z (a f[4, {-1, b2 + 1, b3, b4>] 

- f[4, {-1, b2, b3, b4>]); 

f[4, {bl_, -1, b3_, b4_>] := f [4, {bl, -1, b3, b4>] = 
f [4, {bl + 2, -1, b3, b4>] 

+ z(a“(-l - bl) f[2, {b3 - 1, b4>] 

- f[4, {bl + 1, -1, b3, b4}]); 

f [4, {bl_, b2_, b3_, b4_>] := f[4, {bl, b2, b3, b4>] = 

f [4, {bl + 2, b2, b3, b4>] 

+ z (a~(-bl - 1) f[4, {-1, b2 + 1, b3, b4>] 

- f [4, {bl + 1, b2, b3, b4>]) 

f[5, {-1, - 1, b3_, -1, b5_>] := f [5, {-1, -1, b3, -1, b5}] = 
a~(-b3) f[l, {b5 - 1>] + z (f[3, {b3 - 1, -1, b5>] 

- a~(-l) f[3, {b3, -1, b5>]); 

f [5, {-1, -1, b3_, b4_, b5_>] := f[5, {-1, -1, b3, b4, b5>] = 
a~(-b3) f[3, {-1, b4 + 1, b5>] + z (f[3, {b3 - 1, b4, b5>] 

- a~(-l) f[3, {b3, b4, b5>]) ; 

f[5, {-2, -1, b3_, b4_, b5_}] := f[5, {-2, -1, b3, b4, b5}] = 

a~(-1) f[3, {b3, b4, b5>] + z (a f[3, {b3 - 1, b4, b5}] 

- f[5, {-1, -1, b3, b4, b5>]); 

f[5, {-1, -2, b3_, b4_, b5_>] := f[5, {-1, -2, b3, b4, b5}] = 

f [3, {b3 - 1, b4, b5>] 

+ z (f [5, {-1, -1, b3, b4, b5>] - a~(-2) f[3, {b3, b4, b5>]); 
f [5, {-1, b2_, b3_, b4_, b5_>] := f[5, {-1, b2, b3, b4, b5>] = 

f[5, {-1, b2 + 2, b3, b4, b5>] + z (f[5, {-1, b2 + 1, b3, b4, b5>] 

- a~b2 f[3, {b3, b4, b5}]); 

f[5, {-2, b2_, b3_, b4_, b5_>] := f[5, {-2, b2, b3, b4, b5>] = 

a~b2 f[3, {b3, b4, b5}] + z (a f[5, {-1, b2 + 1, b3, b4, b5}] 

- f[5, {-1, b2, b3, b4, b5>]); 

f[5, {bl_, -1, b3_, b4_, b5_>] := f[5, {bl, -1, b3, b4, b5}] = 
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f[5, {bl + 2, -1, b3, b4, b5>] 

+ z (a~(-1 - bl) f[3, {b3 - 1, b4, b5>] 

- f[5, {bl + 1, -1, b3, b4, b5>]); 

f[5, {bl_, b2_, b3_, b4_, b5_>] := f[5, {bl, b2, b3, b4, b5>] = 

f [5, {bl + 2, b2, b3, b4, b5>] 

+ z (a"(-bl - 1) f[5, {-1, b2 + 1, b3, b4, b5}] 

- f [5, {bl + 1, b2, b3, b4, b5}]) 

f [i_, {-1, -1, b3_, -1, b5_, btail_}] : = 

f[i, {-1, -1, b3, -1, b5, btail}]= a~(-b3) f[i - 4, {b5 - 1, btail}] 

+ z (f[i - 2, {b3 - 1, -1, b5, btail}] 

- a~(-1) f[i - 2, {b3, -1, b5, btail}]); 

f[i_, {-1, -1, b3_, b4_, btail_}] := f[i, {-1, -1, b3, b4, btail}] = 

a~(-b3) f[i - 2, {-1, b4 + 1, btail}] 

+ z (f[i - 2, {b3 - 1, b4, btail}] 

- a~(-l) f[i - 2, {b3, b4, btail}]); 

f[i_, {-2, -1, b3_, btail_}] := f[i, {-2, -1, b3, btail}] = 

a~(-l) f[i - 2, {b3, btail}] + z (a f[i - 2, {b3 - 1, btail}] 

- f[i, {-1, -1, b3, btail}]); 

f[i_, {-1, -2, b3_, btail }] := f[i, {-1, -2, b3, btail}] = 

f[i - 2, {b3 - 1, btail}] + z (f[i, {-1, -1, b3, btail}] 

- a~(-2) f[i - 2, {b3, btail}]); 

f[i_, {-1, b2_, btail_}] := f[i, {-1, b2, btail}] = 

f[i, {-1, b2 + 2, btail}] + z (f[i, {-1, b2 + 1, btail}] 

- a~b2 f[i - 2, {btail}]); 

f[i_, {-2, b2_, btail_}] := f[i, {-2, b2, btail}] = 

a~b2 f[i - 2, {btail}] 

+ z (a f[i, {-1, b2 + 1, btail}] - f[i, {-1, b2, btail}]); 

f[i_, {bl_, -1, b3_, btail_}] := f[i, {bl, -1, b3, btail}] = 

f[i, {bl +2, -1, b3, btail}] 

+ z (a~(-l - bl) f[i - 2, {b3 - 1, btail}] 

- f[i, {bl + 1, -1, b3, btail}]); 

f[i_, {bl_, b2_, btail_}] := f[i, {bl, b2, btail}] = 

f[i, {bl + 2, b2, btail}] + z (a~(-bl - 1) f[i, {-1, b2 + 1, btail}] 

- f[i, {bl + 1, b2, btail}]). 
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